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We present a formalism for analyzing a full-sky temperature and polarization map of the cosmic 
microwave background. Temperature maps are analyzed by expanding over the set of spherical 
harmonics to give multipole moments of the two-point correlation function. Polarization, which 
is described by a second-rank tensor, can be treated analogously by expanding in the appropriate 
tensor spherical harmonics. We provide expressions for the complete set of temperature and po- 
larization multipole moments for scalar and tensor metric perturbations. Four sets of multipole 
moments completely describe isotropic temperature and polarization correlations; for scalar metric 
perturbations one set is identically zero, giving the possibility of a clean determination of the vector 
and tensor contributions. The variance with which the multipole moments can be measured in 
idealized experiments is evaluated, including the effects of detector noise, sky coverage, and beam 
width. Finally, we construct coordinate-independent polarization two-point correlation functions, 
express them in terms of the multipole moments, and derive small-angle limits. 
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I. INTRODUCTION 



With the advent of a new generation of balloon-borne and ground-based experiments ]l| and satellite missions [^,^) , 
the cosmic microwave background (CMB) will provide an unprecedented window to the early Universe. In addition 
to determining the origin of large-scale structure, it has been argued that CMB temperature maps may determine 
cosmological parameters and the ionization history of the Universe, and perhaps probe long-wavelength gravitational 
waves U-|| . 

Any mechanism which produces temperature anisotropies will invariably lead to polarization as well Jici-O]. Tem- 
perature fluctuations are the result of perturbations in the gravitational potentials, which contribute directly to the 
fluctuations via gravitational redshifting (the Sachs- Wolfe effect ) and which drive acoustic oscillations of the pri- 
mordial plasma Q . These processes result in temperature fluctuations which are the same order of magnitude as the 
metric perturbations. In contrast, polarization is not directly generated by metric perturbations: a net polarization 
arises from Compton scattering only when the incident radiation field possesses a non-zero quadrupole moment p6| , |l3[ , 
but only monopolc and dipolc fluctuations are possible as long as the photons in the Universe remain tightly coupled 
to the charged electrons. Polarization is only generated very near the surface of last scattering as the photons begin to 
decouple from the electrons and generate a quadrupole moment through free-streaming p7| . Since by this time most 
of the electrons have recombined into neutral hydrogen, the number of scatterers available to produce polarization is 
reduced, so CMB polarization fluctuations are characteristically at a part in 10 6 , an order of magnitude below the 
temperature fluctuations. 

A polarization map will provide information that complements that from a temperature map. For example, po- 
larization may help distinguish the gravitational-potential and peculiar-velocity contributions to the acoustic peaks 
in the temperature-anisotropy power spectrum [ pd| . In models with reionization, some of the information lost from 
damping of temperature anisotropies will be regained in the polarization spectrum |tq] . Perhaps most importantly, 



* kamion® phy s . Columbia, edu 
t akosowsky @ cf a .harvard .edu 
^stebbinsQfnal.gov 



1 



the density-perturbation and gravitational-wave or vorticity contributions to the anisotropy can be geometrically 



decomposed with a polarization map 16 19|-|21[]. Furthermore, although these non-scalar signals are expected to be 
small, they will not be swamped by cosmic variance from scalar modes (as discussed further below). Detection of 
gravity waves is important for testing inflation and for learning about the inflaton potential which drove inflation [B2| . 

Realistically, detection will present a significant experimental challenge. Current results limit the magnitude of 
linear polarization to roughly a part in 10 5 [^3| . Experiments being planned or built will improve sensitivities by at 
least an order of magnitude . The MAP satellite will make polarized measurements of the entire microwave sky in 
around a million pixels with a precision of around one part in 10 5 per pixel Q . If CMB polarization is not discovered 
by a ground or balloon experiment in the next four years, this satellite will almost certainly make the first detection. 
The COBRAS/SAMBA satellite may also make polarized measurements ||. These experimental prospects, as well 
as the theoretical considerations above, motivate the analysis presented in this paper. 

Previous theoretical treatments of CMB polarization have relied on a small-angle approximation, which is valid 
when considering patches of the sky small enough to be approximated as flat. Upcoming polarization maps will require 
a more sophisticated formalism. In this paper, we develop in detail a description of polarization on the full sky. The 
Stokes parameters conventionally used to describe polarization are not invariant under rotations of the coordinate 
system used to describe them, unlike temperature fluctuations, but rather transform as a second-rank tensor p3| . 
By expressing the polarization in terms of a complete, orthonormal set of tensor basis functions on the celestial 
sphere, power spectra and correlation functions which are independent of the coordinate system can be constructed. 
Earlier work on small patches of the sky chose a particular reference coordinate system which completely defines the 
polarization but obscures the physical interpretation of the polarization pattern. Also, the signal from vector and 
tensor perturbations is expected to contribute to CMB polarization primarily at large angles on the sky through 
gravitational effects, so the correct full-sky analysis is essential. 

Our formalism is stated in terms of differential geometry on the sphere, using a notation widely used in general 
relativity. Similar calculations have recently been performed by Seljak and Zaldarriaga p0|,|21|, using spin-weighted 



spherical harmonics 25 1 . Although the formalisms employed differ substantially and the calculations are quite lengthy, 
we have verified that the end results are equivalent where they overlap, giving us confidence both are correct. 

After a brief review of Stokes parameters, the next Section defines the tensor spherical harmonic basis functions 
and gives useful explicit expressions and formulae for decomposing a polarization map into its harmonic components. 
Section III covers the statistics of the expansion coefficients of the temperature and polarization harmonics, deriva- 
tions of variance estimates for the various multipole moments in idealized experiments, and a recipe for simulating 
a combined polarization and temperature map given theoretical angular power spectra. Section IV derives exact 
expressions for all of the multipole moments from scalar and tensor metric fluctuations, expressed in terms of the 
conventional Fourier components of radiation brightnesses. Section V then treats two-point correlation functions of 
the Stokes parameters in a coordinate-independent manner and expresses the multipole moments and correlation 
functions in terms of each other. We also reproduce flat-sky results by taking small-angle limits and make an explicit 
connection with earlier work in particular fixed coordinate systems. Finally, a summary and discussion section briefly 
considers detection prospects for various polarization signals. A pair of mathematical appendices collect results from 
differential geometry on the sphere and useful identities of Legendre polynomials and spherical harmonics. 

II. DESCRIPTION OF POLARIZATION 

A. Review of Stokes Parameters 

The cosmic microwave background is characterized completely by its temperature and polarization in each direction 
on the sky (assuming its frequency spectrum is a perfect blackbody) . Polarized radiation is described in terms of the 
Stokes parameters Q, U, and V p^j. For a monochromatic electromagnetic wave of frequency ujq propagating in the 
z-direction, the components of the wave's electric field vector at a given point in space can be written as 

E x = a x (t) co$[u} t - x {t)\, E y = a y (t) cos[u t - 6 y (t)]. (2.1) 

If these two components are correlated, then the wave is said to be polarized. The Stokes parameters are defined as 
the time averages 

I = (a 2 x ) + (a 2 y ) , (2.2) 

Q =(<£)-(<%), (2.3) 

U=(2a x a y cos(6 x -9 y )), (2.4) 

V= (2a x a y sm(6 x - 6 V )) . (2.5) 
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The parameter / gives the radiation intensity which is positive definite. The other three parameters can take either 
sign and describe the polarization state. For unpolarized radiation, Q = U = V = 0. The Stokes parameters are 
additive for incoherent superpositions of waves, which makes them natural variables for describing polarized radiative 
transport. 

In most applications polarization is measured in units of intensity; however it is conventional and convenient 
when studying the CMB to express polarization in terms of the difference in brightness temperature of a particular 
polarization state from that of the mean brightness temperature of the CMB. The rationale for this convention comes 
from the well-known result that the spectrum of polarization induced in the CMB is exactly the same as a temperature 
anisotropy, so in brightness temperature units the polarization should be independent of frequency. 

The Stokes parameters / and V describe physical observables and are independent of the choice of coordinate 
system. However, Q and U describe orthogonal m odes of linear polarization and depend on the axes in relation to 
which the linear polarization is defined. From Eqs. ( |2.5| ), it is easy to show that when the coordinate system is rotated 
by an angle a, the same radiation field is now described by the parameters 

Q 1 = Qcos(2a) + t/sin(2a), 

U' = -Qsin(2a) + t/cos(2a). (2.6) 

Stated another way, under rotations of the coordinate system around the direction of propagation, the Q and U Stokes 
parameters transform like the independent components of a two-dimensional, second rank symmetric trace-free (STF) 
tensor. Thus we can equally well describe the linear polarization state by a polarization tensor Tab, which coincides 
with the photon density matrix fl3]j . 



B. Scalar and Tensor Harmonic Expansions 

Suppose we have an all-sky map of the CMB temperature T(n) and polarization tensor (n). The polarization 
tensor is a 2 x 2 symmetric (V a b = Vba) and trace-free (g ab V a b = 0) tensor, so it is specified by two real quantities. 
Given the Stokes parameters Q and U measured in any coordinate system, we can construct V a b- For example, in 
spherical polar coordinates, (9,4>), the metric is g a b — diag(l, sin 0) and 

1 / Q(n) -U(n)sm6 \ 

The factors of sin# must be included since the coordinate basis for (9,<fi) is an orthogonal, but not an orthonormal 
basis. (For more details of differential geometry on the 2-sphere, see Appendix A.) The Compton scattering process 
which thermalizes the CMB and generates polarization cannot produce any net circular polarization p7|] ; thus we 
expect V — for the microwave background and do not consider the V Stokes parameter further. Note the spherical 
polar coordinate system adopted in this paper gives an outward direction for the z-axis, which is opposite the radiation 
propagation direction. The convention with the z-axis in the direction of propagation is sometimes used, particularly 
in Ref. fisjj ; this leads to the opposite sign for the U Stokes parameter, but all results are unchanged. 

In the usual way, we can expand the temperature pattern T(n) in a set of complete orthonormal basis functions, 
the spherical harmonics: 



T(n) 
To 

where 



T 

a (lrn) 



1 + EE 4m)Y(lm)(&) (2.8) 



l—l m— - 



i- J d&T{n)Y$ m) (n) (2.9) 



are the temperature multipole coefficients and To is the mean CMB temperature. The / = 1 term in Eq. is 
indistinguishable from the kinematic dipole and is normally ignored. 

Similarly, we can expand the polarization tensor in terms of a complete set of orthonormal basis functions for 
symmetric trace-free 2x2 tensors on the 2-sphere, 



Vab(n) 



OO t 

EE [ fl ° m ) r (L)a fc (fi) + ^ m )^(L )a6 (n)] , (2.10) 
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where the expansion coefficients are given by 



4m) =yJ d ^a b (n)Y ( f^ b *(h), af lm) = ±-J dnV ab (n)Y™*(n), (2.11) 

which follow from the orthonormality properties 

/ dnY^ )ab {n) F ( «$(n) = J dhY^ )ab {n) = 8 w 8 mm> , (2.12) 



dnY& )ab (n)Y$° b ) (n) = 0. (2.13) 

Note that unlike scalar harmonics, the tensor harmonics only exist for / > 2 [f28f . 

The basis functions Yff m <. ab (n) and Yfi m <. ab (ii) are given in terms of covariant derivatives of the spherical harmonics 

by II 

Y (f m ) ab = N[ (Y(l m ):ab ~ ^ 9abY(l m ).J\ , (2.14) 



and 



^Y(lm):ac eC b + Y^ lm y bc e C a ^ , 



Ya m )ab — "7T ( Yn m y ac e c b + Y, lm y bc £ c a ) , ( 2 -15) 



(lm)ab <2 

where e a b is the completely antisymmetric tensor, the ":" denotes covariant differentiation on the 2-sphere, and 



is a normalization factor. 

The existence of two sets of basis functions, labeled here by "G" and "C" , is due to the fact that an STF 2x2 
tensor is specified by two independent parameters. In two dimensions, any STF tensor can be uniquely decomposed 
into a part of the form A- ab — (l/2)g ab A :c c and another part of the form B :ac e c b + B :bc e c a where A and B are two 
scalar functions. This decomposition is quite similar to the decomposition of a vector field into a part which is the 
gradient of a scalar field and a part which is the curl of a vector field; hence we use the notation G for "gradient" and 
C for "curl". Since the Y(; m )'s provide a complete basis for scalar functions on the sphere, the Y^f m ^ ab and Y9\ ab 
tensors provide a complete basis for G-type and C-type STF tensors, respectively. This G/C decomposition is also 
known as the scalar /pseudo-scalar decomposition p8| |. 

Incidentally, these tensor spherical harmonics are identical to those which appear in the theory of gravitational 
radiation [§9],f30| • The propagating degrees of freedom of gravitational field perturbations are described by a spin-2 
tensor. Computing the flux of gravitational radiation from a source requires the components of the gravitational field 
tangent to a sphere around the source which are induced by the motions of that source. Our G harmonics are often 
— but not always |3(| — referred to as having "electric-type" parity, since an electric field can be written as the 
gradient of a scalar. Likewise, our C harmonics have "magnetic-type" parity since they are the curl of a vector field. 
The two varieties of harmonics also corre spon d to electric and magnetic multipole radiation. 



Integration by parts transforms Eqs. (2.11) into integrals over scalar spherical harmonics and derivatives of the 
polarization tensor: 



'!/'„,, = Y I dnY { * lm) (h)V ab ab (n), (2.17! 



c = W 
To 



a {lm) 



J dnY { * lm) (n)V a b ac (n)e c b , (2.18) 



where the second equation uses the fact that e :c = 0. These forms are useful for theoretical calculations of the 
multipole moments. We don't recommend taking second derivatives of real data! Since T and V a b are real, all of the 
multipole must obey the reality condition 

= (-l) m 4— )• ( 2 - 19 ) 

where X = {T,G,C}. 
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C. Explicit Form of the Harmonics 



In spherical polar coordinates, (9, </>), the tensor spherical harmonics are given explicitly by p^j30[ |, 

q ,^_Nif w (irn){n) X {lm) (h)sm9 \ 

Y (l m )ab\ n ) — TT [ v ( " \ • ft ii, i-\ ■ 2ft (2.20) 
2 \X {lm) (n)sm9 -W (lm) (n)sm 2 9 J 



and 



= 4 1 w . a v ~"'\: . 2 j > ( 2 - 21 ) 

2 ^ w (Im) ( n ) sin e X (lm) (n)sin 2 0/ 



iV ; / - x (lm)(&) ^(im)(n)sin( 



where 



and 



/ d 2 d m 2 \ ( d 2 \ 

«W fi ) = y W ~ cot ^ + —2- e ) Y dm)(n) = 2^— -1(1 + 1) j r (M (n), 



2 £777. / ^ \ 



(2.22) 



(Note that this definition of X^ m j(h) differs from that in Ref. |30| by a factor of sinfl.) The exchange symmetry 
{Q, U} <-> {U, —Q} as G^C indicates that Y(f m } ab and Yn m ) a ), represent polarizations rotated by 45°. By evaluating 
the derivatives, these functions can be written 



W( lm) (*) = - 2 J 2 4t i [s + rji G ^) (c ° sg)e " n "' (2 - 24) 



-\ t, h G/j i cost) e 

y 4tt (l + m)\ ( ,m > 1 > 



iX (lm) (n) = -2J - ), ; j; <Wcosfl)e' m » (2-25) 



where the real functions G^ m ) are defined by j25 



G+ (cos0) = - (K^r + 1)) + (J + r^^-P^cosQ), (2.26) 

1 ; V sm 9 2 J sin 



G^ m) (cos 0) = ((i - 1) cos 0if (cos 9) -{I + m)P£ 1 {x)). (2.27) 

These expressions will be useful for the correlation functions in Section V, and for simulating maps and data analysis. 

In linear theory, scalar perturbations can produce only G-type polarization and not C-type polarization. On the 
other hand, tensor or vector metric perturbations will produce a mixture of both types . Heuristically, this 

is because scalar perturbations have no handedness so they cannot produce any "curl", whereas vector and tensor 
perturbations do have a handedness and therefore can. Observation of a non-zero primordial component of C-type 
polarization (a nonzero a^ m )) in the CMB would provide compelling evidence for significant contribution of either 
vector or tensor perturbations at the time of last scattering. 

Given a polarization map of even a small part of the sky, one could in principle test for vector or tensor contribution 
by computing the combination of derivatives of the polarization field given by "P ab : b c e c a which will be non-zero only 
for C-type polarization. Of course, taking derivatives of noisy data is problematic. We discuss more robust probes of 
this signal below. 
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III. STATISTICS OF THE MULTIPOLE COEFFICIENTS 



A. Statistical Independence of the Coefficients 

We now have three sets of multipole moments, aji m y a< n m \i an d ®?i m y which fully describe the tempera- 
ture/polarization map of the sky. Statistical isotropy implies that 

/ T * T \ y^T c r / G * G \ r r 

\ a (lm) a (l'rn') ) = °W°mm>, \ a (lm) a (l'm') / = °U'°mm>, 

C * C \ /^«C r s" / X * G \ z^fTG r .c 

a (lm) a (l'm') ) — Oll'Omm'i \ a (lm) a (l>7n') / = °U'"mm', 

( a Jm) a (i'm')) — Q TC ^«"^""n'; (^ a flm) a fl'm') ) = Cf C 8w 5 mm > , (3.1) 

where the angle brackets are an average over all realizations. For Gaussian theories, the statistical properties of a 
temperature/polarization map are specified fully by these six sets of multipole moments. In fact, the scalar spherical 
harmonics Y(; m ) and the G tensor harmonics Yn m \ a ^ have parity (—1)', but the C harmonics Y^ m ^ ab have parity 

(— Therefore, symmetry under parity transformations requires that C; TC = C GC = 0, which will also be 
demonstrated explicitly in the following Section. Measurement of nonzero cosmological values for these moments 
would be quite extraordinary, demonstrating a handedness to primordial perturbations. In practice, these two sets of 
moments can be used to monitor foreground emission. Furthermore, as mentioned above and demonstrated explicitly 
in Sec. IV, Cp = for scalar metric perturbations p9| , p0| . At small angular scales where the contribution from 
tensor and vector perturbations is expected to be negligible, Cf can also be pressed into duty as a foreground 
monitor. Exact expressions for these multipole moments in terms of the photon brightnesses usually calculated by 
early-Universe Boltzmann codes are derived below. 



B. Map Simulation 



For the case of Gaussian statistics, realizations of temperature/polarizations maps are easy to generate using 
standard techniques. Since the only cross-correlation between mode coefficients, given by C ; TG , correlates only aj lm ^ 
and a G m ) with the same / and m, the total correlation matrix is block diagonal with the largest blocks being only 
2x2 matrices. In particular, set 



\lm) 



Ci (Q T ) 1/2 . 



Ci 



(c T } l/2 
a flm) = C3 (Cp) 



1/2 



(3.2) 



where for each value of I and m > choose three complex numbers, (£i, £2, C3) ; drawn from a Gaussian distribution 
with unit variance, i.e., both v / 2Re(Ci) and V21m(£i) are drawn from a normal distribution. For m — the same 
equations hold but the Ci should be real and normally distributed; for m < the coefficients are given by Eq. ( |2. 19 ). 
Note that in all cases CpCf > (C Z TG ) . This set of coefficients can be combined with Eqs. ($7$, ( plcj ), ( |2~23| ), and 
( [2.22 ) to obtain the explicit expressions 



Q(n) = 27Vn) = T ]T £ iV ; [a G m) IU (im) (n) - af lm) X (lm) (n)] 

1=2 m=-l 

00 I 

U(h) = -2csc07V(n) =-T J2Yl ^["(M^Wn) + af lm) W {lm) (h)] 



(3.3) 



with Wn m \ and Xn m \ given by Eqs. (2.24) and (2.21:). Note that polarization maps are traditionally plotted as 



headless vectors with amplitude (Q 2 + C/ 2 ) 1 / 2 and orientation angle (1/2) tan 1 (U/Q). 
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C. Estimators 



One of the the main uses of a temperature/polarization map will be to determine the multipolc moments with the 
best possible accuracy. From a full-sky CMB temperature map, we can construct the following rotationally-invariant 
estimators (denoted by a ~) for the multipolc coefficients: 



c?= E 



\a 



{It, 



I n.9. 



21 + 1' 

2 



c?= E 



(lm) 



m——l 



21 + 1' 



l 

E ■ 

■m— — l 


a {lm) 


2 


21+1' 


I 

= E 

rn — — 


T * 
07, > 


a?, n 


21 


+ 1 



(3.4) 



Note that Eq. ( p. 19 ) guarantees that the Cj G will be real. 

When averaged over the sky (denoted by an overbar), the mean square temperature anisotropy after subtracting 
the dipole is 



1=2 



4tt 



l • 



(3.5) 



and the mean square polarization is 



V 2 = Q 2 + U 2 = 2V ah Vab = V G + Vc 



1=2 



1=2 



(3.6) 
(3.7) 



Even if no single CP or CP gives a significant signal, combining different l's as in Pq 2 or Vq 2 can give a statistically 
significant signal. 



D. Cosmic and Pixel-Noise Variance 

The averages in Eqs. ( |3.l| ) are over an ensemble of universes drawn from a theoretically defined statistical distri- 
bution, or assuming ergodicity, a spatial average over all observer positions in the Universe. However, we can only 
observe a single realization of the ensemble from a single location. Therefore, even if we had an ideal (full-sky coverage, 
no foreground contamination, infi nite angular resolution, and no instrumental noise) experiment, the accuracy with 
which the estimators in Eqs. (^J) could recover the multipole moments would be limited by a sample variance known 
as "cosmic variance." Furthermore, a realistic experiment may have limited sky coverage and angular resolution and 
some instrumental noise. In this Section, we calculate the cosmic variance with which the multipolc moments can be 
recovered. We also calculate the variance due to finite sky coverage, angular resolution, and instrumental noise in an 
idealized experiment. To do so, we adopt a simplified model in which we assume a pixelized map in which the noise 
in each pixel is independent and Gaussian distributed after foregrounds have been successfully subtracted. In many 
respects, our derivation follows that in Ref. |32|, and our results agree with those in Ref. 

We must first determine the contribution of pixel noise to each multipole moment, and we begin with the temperature 
moments. Consider a temperature map of the full sky, T map (n), which is pixelized with N p - iK pixels. If we assume 
that each pixel subtends the same area on the sky then we can construct multipole coefficients of the temperature 
map using 

,t f , /T ma P(n)\ , , 1 ^ 4tt man , N , . 

dji m ) =jdn ^ > ) y (Im) (fi) ~ - J2 T T Y(im)(Zj), (3.8) 

where T™ ap is the measured temperature perturbation in pixel j, and fL, is its direction. The difference between d(; m ) 
and <Z(i m ) is that the former includes the effects of finite beam size and detector noise, i.e. dn m \ is measured. The 
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extent to which the approximate equality fails is the pixelization noise, which is small on angular scales much larger 
than the pixels. In many of the proposed experiments, which oversample the sky when compared to their beam, 
there will be no loss of information by ignoring scales close to the pixel scale. The observed temperature is due to a 
cosmological signal and a pixel noise, Tj nap = Tj + T^ olse . If we assume that each pixel has the same rms noise, and 
that the noise in each pixel is uncorrelated with that in any other pixel, and is uncorrelated with the cosmological 
signal, i.e. (T^T? ™) = T 2 (o-J ix ) 2 <% and (T t T™ isc ) = 0, then 

/ jT jT* \ / T T* \ , / T, noise / T,noise\ * \ 
(<V) C W)/ = \ a (lm) a (l'm')/ + { a (im) {"(I'm') ) ) 

= \Wh 2 C?6 w 5 mm , + (agf ° (a^ c )*), (3.9) 

where we have written the expectation value of the cosmological signal in terms of that predicted by theory, C ; T , 
multiplied by IW^I 2 which accounts for beam smearing. Typically the beam is approximately Gaussian in shape, 
corresponding to the window function w exp(— Z 2 ct 2 /2), where er b = 0f w h m /V81n2 = 0.00742 (0f w hm/l°) and 
^fwhm gives the full- width at half maximum. 
The second term in Eq. (3.9) is 

1 ^V4^ 2 



<«S) iS0 ( a S ; ™t) ) = j* E E (wr) (^ oise Tr ise >i(M(fii)^ m o( fi i) 

i—l j — 1 ^ P 1X / 

= fe) 2 E j^ y (im)(^0fym')( fi i) 



> 2 



47T (crj iv , 

y ' S w S rnm ,. (3.10) 



Therefore, the moments measured by the map are distributed with a variance |32[ 

>2\ 



<$m)4m')) = \Cl\W?\ 2 + j S W S mm , (3.11) 

Now we move on to the noise contribution to the polarization moments. We will assume the instrumental noise 
in the polarization measurements is isotropic, the same for all pixels, uncorrelated with the noise in the anisotropy, 
which in terms of the Stokes parameters requires 

(Qf oisc Qf lsc ) = (U™ isc U™ isc ) = :T Vpix) 2 % (Qr isc U™ isc ) = (Q? oisc T™ isc ) = (u™ isc T™ isc ) = 0. (3.12) 

We denote the polarization tensor describing the noise in pixel i by V^ lsc (iii) . The previous equations are equivalent 
to the coordinate independent equation (see ref p8|) 

(^f sc (fl,) Kf c (fi,)) = \ tf{^) 2 {g ac g b<i - e oc e M ) <% (Prffi.) T"° isc (n,)) = 0, (3.13) 



The mode coefficients for the noise, defined as in Eq. (2.11), will have a correlation matrix 



X, noise * X' , noise * \ / 
a {lm) a (l'r 



E E CPai{&i)Vcd(&i)) X, X' G {G, C} (3.14) 

»=1 3=1 



so using Eqs. (flf), ( ^13|) , and @ we find 

C. P \ 2 Npix 

lVp~ J E 21 W) ( n i) Y (/'m')«^ ni ) = iVp-Kix) S w S mm >8xx> (3.15) 



and of course 

'\lrn) Ll (l'm') / ~ V^im) "(I'm') 



G. noise* T,noisc\ / C, noise* T,noisc\ n / i n\ 
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Thus instrumental noise contributes equally to the measured G- and C- polarization components but introduces no 
cross-correlation between them. 

Collecting the results, if d x m -, (for X = {T, G,C}) are the multipole coefficients for the map (signal plus noise), 
their variances will be 

(4™)4' m <)) = (Kl 2 Q XX ' + w£,)8u>S mm > = Df x '5 lv 5 mm ,, (3.17) 

where 

-1 _ 4?r ( a Jix) 2 -1 -1 _ -1 _ 2?r Kix) 2 ,_ 1H » 

w T = — — E , Wgg = w cc = Wp = — , (3.18) 

and w^^, = for X ^= X'. The quantities w _1 are inverse statistical weights per unit solid angle, a measure of 
experimental sensitivity independent of pixel size [Q. (Note that our Wp 1 differs by a factor of two from that in Ref. 
pjj which is consistent with our Cp and Cf differing from their E and B moments by a factor of two.) 
Estimators for the multipole moments CMB power spectra, C xx , are 

Cr' = (Dr'-w^,)\Wt\- 2 , (3.19) 



where 

l d x* d x' 

DfX , = ^ JgUjpL. (3.20) 

m— — / 

Since the estimators for each D xx (and therefore for C xx ) are constructed from only 21 + 1 multipole coefficients, 
each C xx can be recovered only with a finite sampling variance, known as cosmic variance. In addition, the six 
different sets of measured moments are constructed from the three sets of dn m ) coefficients, leading to some covariance 
between the moments. These variances can be described with a (6 x 6) covariance matrix, 



■=AA' 



^DfDf)-D?Df)\Wl>\- 4 (3-21) 
for A= XX'. 

We now calculate the entries of this matrix. Recall that if Xi are Gaussian random variables with variances 
(xiXj) = crjj, then {x^xf) = crf^ + 2cr?-, and (x^Xj) = Sa^crfj. For X={T,G,C}, 



( D fX Df X'\ = £ 



(21 + iy 



E 



; (2i + 1) 2 

XX nXX' 



D XX D X ' x '(l - 5 mm ,) + ( D XX D X ' X ' + 2 (D xx '^ 



Dr + {2l + l)2 (d XX ') . (3.22) 
The diagonal elements for TT, GG, and CC are thus 

Sxxxx = ^ (C xx + w^\Wr\- 2 f , (3.23) 

and the off-diagonal elements are 



™ <: 2F—\" 

Stt,cc = 0, 

Sgg,cc = 0. (3.24) 
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For the diagonal TG component, 



E 

mm' 

(2T 



^(lm) ^ \lm)^ \lm' ) ^ \lm' ) 
(2/ + 1)2 



v ^po + 2(A TG ) a , V (1 - W) (A™)' 
^ (2Z + l) 2 L> 



(21 + If 



1 



2Z + 1 



(A ) 



Therefore, 



1 



'TG,TG 



2/ + 1 



(3.25) 



(3.26) 



The diagonal covariance-matrix elements for TC and TG are obtained similarly and are nonzero even though C ; TC = 
C GC = 0. Given a map, it should be checked for consistency that there is no statistically significant parity violation. 
Finally, for the off-diagonal TG-TT component, 

(^T G ? T ) = E [(\ d Vm)\ 3 d? lm) )5 mm * + (\dj lm) \ 2 (dj lml) yd? lm/) )(l - 5 mm ,) 



21 + 1 



1 D,D 



TT r,TG 



Therefore, 



and similarly 



"TT,TG 



21 



T [Q TG (Q T +%Wr 2 )] 



(3.27) 



(3.28) 



"GG.TG 



21 



(3.29) 



and this completes the calculation of all nonzero elements of the covariance matrix. 

To determine the precision with which a temperature/polarization map can recover some cosmological parameters, 
e.g., s = {f^oi A, h, fib, •••}> w e can evaluate the curvature, or Fisher information, matrix [ |33"| |, which can be generalized 
from that in Ref. jl]] to 



l A, A' 



AA 



(3.30) 



for A=TT,GG,CC,TG where [S _1 ]^» are elements of the inverse of S. The standard error with which a given 
parameter s, can be recovered after marginalizing over all other parameters is given by the square root of the diagonal 
element i of [a] -1 , assuming a linear dependence of C, on all of the parameters. 

The w^ 1 factors on the right-hand side of Eqs. (3.23), (3.26), and ( 3. 29] ) are those due to instrumental noise. 
However, note that even in an ideal experiment wit h cr P i X = 0, the right-hand sides would still be nonzero, and this 
is the cosmic variance. Eqs. (3.23), (3.26), and (3.29) are valid for a map with full-sky coverage. Realistically, 
however, only a fraction / s ky of the sky will be surveyed, or if the entire sky is surveyed, only a fraction will be used 
in the analysis. Therefore, the the accuracy with which cosmological parameters can be recovered will be degraded 
accordingly. Strictly speaking, harmonic analysis on a cut sky will have to be performed, and this will introduce 
correlations in the errors of multipole moments of different Vs. However, if the entire sky is surveyed, but only a 
fraction of the sky is used for the analysis (e.g., if the Galactic plane has been cut out), then the effect of partial sky 
coverage can be approximated by multiplying the curvature matrix by / s k y (which will increase the standard errors 

— 1/2 

in cosmological parameters by / skv )• If only a fraction of the sky is surveyed, the curvature matrix should still be 
multiplied by / s k y , but note that iVpi x is the number of pixels actually in the map. 
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E. Pixel Noise for a Polarization Map 



How are the temperature and polarization pixel noises related? If the two linear polarization states are always 
given equal integration times, the total number of photons available for the temperature measurement will be twice 
the number available for either polarization measurement. Therefore, 




(3.31) 



However, a crucial difference in overall sensitivity exists between the two current receiver technologies. Coherent 
receivers (i.e. HEMT amplifiers) measure the incoming electric field, while incoherent receivers (i.e. bolometers) 
measure only the incoming total power. In the former case, the signal can be split into two orthogonal linear 
polarizations and the phase information can be retained throughout the entire signal path. If properly designed, such 
a system of receivers can measure the linear polarizations without compromising the temperature measurement, so 
the temperature sensitivity will be determined solely by the amplifier characteristics This is the design strategy 
for the MAP satellite. 

On the other hand, bolometers measure only the instantaneous total power received and do not retain any phase 
information; a filter must be placed in front of the detector for each linear polarization state, discarding half of 
the incoming photons. Thus the temperature sensitivity for a polarized bolometer measurement is only half of that 
for an unpolarized measurement, which collects twice as many photons in the same amount of integration time. 
A compensating factor is that bolometers offer much greater raw sensitivity than HEMT amplifiers. An important 
question facing future bolometer experiments is whether to sacrifice half of the temperature sensitivity for polarization 
information j34|] . If the goal of an experiment is to measure angular power spectra and the temperature measurements 
are dominated by cosmic variance, then polarizing the measurement is obviously advantageous. The answer is as yet 
unclear in cases where cosmic variance is not the controlling factor in the temperature measurement. 



IV. CALCULATION OF THE MOMENTS 



In this Section, we calculate the set of multipole moments defined by Eqs. (3.1), for scalar and tensor metric 
perturbations. Vector metric perturbations make a negligible contribution for inflationary theories, although they are 
generic in defect models; they will be covered elsewhere. The intensity and linear-polarization state of the CMB in 
any given direction is specified by three quantities (the temperature T and the Stokes parameters Q, and {/), giving 
six possible sets of multipole moments, C Z T , Cp, Cp, C ; TG , C ; TC , and Cp c , but as argued above, parity demands that 

C TC = C GC = 

One way to calculate the moments is to rewrite the radiative transfer equations in terms of tensor harmonics |3l| . 
The contribution of each Fourier mode to each multipole moment is then obtained by evolving numerically the coupled 
Einstein and Boltzmann equations for the multipole coefficients o^ [m y Integrating over all Fourier modes then gives 
the multipole moments. This approach has the advantage of being similar in form to the usual moment hierarchy 
formulation of the problem, while keeping the independent modes separated throughout the calculation, giving simple 
power spectrum expressions. A second approach offering computational advantages has been presented in Ref. 
which uses the Stokes parameter evolution equations to write an integral equation solution for the multipole moments. 

Another possibility, which makes contact with previous work on CMB temperature anisotropics and polarization, 
is to express the multipole moments in terms of the usual perturbations to the photon brightness and polarization, 
An and Aqi (see Ref. jl3| for definitions and descriptions), obtained from most current numerical calculations [j35| , 
and this is the approach adopted here. In this Section, we calculate the contributions to these moments from scalar 
and tensor perturbations. With these results, it is straightforward to modify existing numerical codes to obtain all of 
the multipole moments. 



A. Scalar Metric Perturbations 



The simplest calculation of the multipole moments in terms of the photon brightnesses uses the fact that, due to 
statistical isotropy, the contribution of a given k mode to the moments C; depends on its magnitude only and not 
its direction. Therefore, we will consider the contribution of a single k mode, with the coordinate system always 
chosen with z in the k direction, and then integrate over all k at the end. All temperatures, Stokes parameters, 
polarization tensors, and expansion coefficients are functions of k, although we sometimes drop explicitly reference to 
k for notational simplicity. Functions of the real-space coordinate x do not appear in this paper. 
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1. Temperature moments 



We begin with the familiar temperature moments. The temperature anisotropy induced in the direction n on the 
sky by a single k scalar mode is given by Eq. (7.1) in Ref. |l3f| , 



T(k,n) 1 



J2(2j + l)P 3 -(k • n) A^(fc) = 1 + £ E a Jm)( k ) y (ta)(fi)- (4-1) 



T 4 

u j=l 1=2 rn=~l 



The superscripts T and s indicate that we are dealing here with temperature moments from scalar perturbations. 
The A s n (k) are Legendre coefficients of the photon intensity distribution function for scalar metric perturbations [ p^[ . 
The expansion coefficients are given by the inverse transform, Eq. ( |2.9| ), 



T,s 
a (lm) 



= zE( 2 J + V^'iW / ^p J -(k-fi)y ( i ro) (n) 

j=o 

= 7 rA s /i (fc)y ( t m) (k) 
1 



= -y/(2l + l)nA s n (k)S m0) (4.2) 



where we have used Eq. (B12) and taking k = z in the final line. The contribution to C^' s from this mode is then 

m— — I 

The total multipole moment from scalar modes, ^ T ' scalar ^ i s given by integrating over all k: 

-j Wf Cl w_ 8^ 

which is the usual result. 



^scalar = f c T, S(k)= 1 , ,, (lk |A}/(/ . )| ^ !4 .4) 



3. Polarization moments 



Now we move on to the polarization moments produced by scalar perturbations. First we describe the "standard" 
representation of polarization, which is what has usually been computed by Boltzmann codes. In the spherical polar 
coordinates, (9, 4>), the Stokes parameters induced in direction n on the sky by a single k scalar mode in the z direction 
can be obtained from Eq. (7.4) in Ref. [jl3| (the £' there can be chosen zero for k = z) and are 

= I f>' + l ) P & ■ fi ) A Q#) U% n) = 0, (4.5) 

where the Ag ; (fc) are Legendre coefficients of the photon polarization distribution function for scalar metric pertur- 
bations. The polarization tensor at a point n on the sky induced by this scalar mode is thus 

3=0 



^ oo 

= - E( 2 J + l) A Qj(k) M U) ab(k A), (4.6) 



8 ^o 



where the second line defines the tensor ^ i.e., the tensor M?K (with raised indices) takes the form 



M$(k = Z i n)=P j (co S 0)l I _^ 2a ). (4.7) 
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This j-expansion is not an expansion in tensor spherical harmonics and, as we shall see, the harmonic content of the 
M^j tensor, while peaked around I — j, has significant contributions from I far from j. 

We now proceed to re-express the above rep rese ntation of polari zatio n in terms of tensor spherical harmonics. The 
G multipole coefficients of the pattern in Eq. (|4.6|), given by Eq. (2.17), are 



(4.8) 



ag„ s (k) = N t J dnY ( * lm) (n)V ab .. ab (k,n) 

= ^(2j + l)A s Qj (fc) f dhY ( * m) (n)M? b .. ab (h) 



We may compute M" - ab in spherical polar coordinates by substituting Eq. (4/7) into Eq. (A13), obtaining 



M$ :ab = (1 - x 2 ) pV{x) - AxP'^x) ~ 2P j (x), 



= k - fi. 



(4.9) 



This can be simplified using the definition of the associated Legendre polynomials P[ n (x) and the recursion relation 
Eq. (p5|), giving 



-Pf(x)-2(.f+3 + l)Pj{x). 



Pj(x) can be represented as a finite series in Pi{x) using the integral Eq. ( |B10| ), giving finally 



1 



M$ :ab 



■£(l + l)(l + 2)cyfi(£-fi), 



Clj 



2(21 + 1} 



1=3 

I — j even and I < j 



(J + l)(Z + 2) 
I — j odd or I > j 



(4.10) 



(4.11) 



Making use of Eq. (B12) and orthonormality of the spherical harmonics, 

J dhY { l n) (h)M^.. ab (h) = -(I + l)(l + 2 )c ljy J-^- [ 6 m0 , 



and 



where 



A s Gi (fc) 



2(l + l)(l + 2) \ 1/2 ^(2j + l\ 

1(1-1) ) ^ {uTi)^ a q^- 



(4.12) 



(4.13) 



(4.14) 



This infinite sum is, as shown below, equivalent to the finite sum of Eq. (4.24). The contribution to C ; G from this k 
mode is 

^ G ' S(k) = 2TTT £ l^)^ 2 = i^l A ^( fc )l 2 - ( 4 - 15 ) 

TCI — — I 

The moments C/ are rotationally invariant, so assuming statistical isotropy, which guarantees that the different k 

ion c 

d 3 k 



modes are uncorrelated, the total contribution of all scalar modes to C, is 



a 



G, scalar 



(4.16) 



The calculation of Cf ' s is similar with the replacements V ab -. ab -> V ab :a C e c b and Mp yab -> M? b ):ac e° b in Eq. (|48|) 



0) : 



'(i) : 



[cf., Eqs. J2.17|) and (|2.1g|)]. However substituting Eq. (^LTJ) into Eq. ( |A14[ ) yields M ( { lb . ac e c b = since M$ is diagonal 



'00 



and independent of This is just what we expect: Since M?k is even under parity while e c b is odd, the product must 



integrate to zero. Thus for scalar perturbations, a ,s (k)(; m ) = and 
as argued above. 



^C. scalar ^ 

°z — U > 



(4.17) 
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FIG. 1. The basis tensors My'(fc) traditionally used for polarization are discontinuous in the orientation of the polarization 
in the directions on the sky both parallel and anti-parallel to k. The discontinuity is either as depicted in the left panel or as 
depicted in the right panel. Note that we have just switched the sign of the polarization between the two panels. 



3. Polarization in the ±k directions 



As seen from Eq. (4.11), when I is large the coefficients cij at I = j are much larger than "nearby" coefficients, say 
when I = j — 2. If the I = j term dominates, then in the small angle approximation, i.e., for I 3> 1, the approximation 
^Gi(^) ~ Ag ; (fc) is valid. However, this is a priori not a very good approximation since the contribution from 
the terms with I + 2 < j < 21 to the sum of Eq. ( 4.14 ) comes to nearly as much as the contribution from the j = I 
term (although some cancellation may result from sign changes in Aq ; ). The explanation is that it takes the sum of a 

large number M^(n) to represent Y^ b (n). This behavior is expected for the simple reason that while the Y^^'s are 
smooth functions, the M^'s are not: for n = ±k, i.e., when 6 — and it (which are singular points of the spherical 
polar coordinate system), does not go to zero (since Pj(l) ^ and Pj(— 1) 7^ 0). Instead the amplitude of the 
polarization approaches a constant but its direction varies discontinuously as illustrated in Fig. |l|. To represent this 
discontinuous behavior as a sup erposition of smooth functions requires a large number of terms. In fact, the only 
reason why the sum of Eq. (4.11) does not contain an infinite number of terms is because it includes either I = if j 
is even or I = 1 if j is odd, neither of which are part of the basis of harmonic STF tensors. 

While the basis functions are discontinuous, the physical polarization pattern must be continuous. The 

polarization will in fact be zero in the directions n = ±k as can be seen directly from the form of the Boltzmann 
equations [T^ ] . Thus the solution of the Boltzmann equation will obey 



V ab (k,+k) cx ^(2j + l)A s Qj (k) = V ab (k, -k) ex J2(-l) j (2j + l)A s Qj (k) = o 
3=0 



(4.18) 



3=0 



which in turn implies 



odd 



]T(2j + l)A^-(fc) = £(2j + l)A s Qj (k) = 0. 

j>0 j>l 

Assuming a continuous polarization pattern and substituting Eqs. ( |4.6| ) and (4.11) into Eq. ( 4.12| ) gives 

/oo even 
dnrf 00) (n)PW(k,n) ex £(2j + 1) C <Wq#) = + ^ A Q#) = °> 



(4.19) 



3=0 



j>0 
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/ 



odd 



dhY{ w) (h)V ab :ab (k,n) « £(2j + l) Clj A s Qj (k) = £(2j + 1) A° Qj (k) = 0. 

J=0 j>l 



Thus the decomposition of P ab :a b actually has no / = or / = 1 content. 
Define the coefficients 



I 2 -I 



bu = { 



(/ + l)(/ + 2) 
2(21 + 1) 
+ + 2) 



I 



/ — j even and < j < I 
I — j odd or j > I 



which have the property 
(/ + !)(/ + 2) 



2(2/+ 1) 



hj (2j + 1) A s Qj (fc) + J2 cy (2J + 1) A s Qj (fc) 

J=0 j=J 



^(2j + l)A s Qj (fc) /even, 



j>0 
odd 



^(2j + l)A s Qj (fc) /odd. 

3>1 



Since these sums are zero for smooth (cosmological) polarization patterns, we may use the equality 



(4.20) 



(4.21) 



(4.22) 



£cy(2i + l)A 8 Qi (fe) 

3=1 



]TM2j + l)A^(fc) 

3=0 



in Eq. ( pl| ) to obtain 



A b Gi (fc) = 



2(/ + !)(/ + 2) 
Z(J-l) 



1/2 ! 

E 



2j_+l 
21 



r ) fc«A" 0i (Jfe). 



(4.23) 



(4.24) 



We have transformed an infinite sum into a finite sum. While these finite sums are still somewhat cumbersome, 
they are significantly less complicated than previous expressions for moments of Stokes parameters obtained in the 
small-angle limit [[l3 14 1. 

The fact that continuity demands zero polarization at n — ±k is reflected in that both Y^ ab and Y^ ab are zero 
at 9 — and n, or equivalently that G^(±l) = 0; continuity in fact requires G^ m ^(±l) = for all m except when 

m = ±2. Since the m = ±2 terms correspond to tensor perturbations when k and z are aligned, the polarization 
from tensor perturbations does not give large sums, as evidenced below. However, a similar treatment of vector 
perturbations (m = ±1) leads to the same sort of infinite sums encountered for scalar perturbations. 



4- Cross moments 



From Eqs. (4.2) and (4.13), the contribution to the TG cross moments from a single scalar mode in the k = z 
direction is 



c ' TG ' S(k) = 27+T E ( a £)( k ))* a S( k )> 



(4.25) 



so integrating over all k gives 



C, 



TG, scalar 
I 



J (2^ TG ' S(k) = lib / ^dkAUk^ik). 



(4.26) 



The vanishing of a^-j(k) also demonstrates explicitly that the moments C; TC ' S = C ; GC ' S = 0, as argued above. 
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B. Tensor Metric Perturbations 



1. Temperature moments 



The calculation of tensor moments proceeds in an analogous fashion. Recall, however, that tensor perturbations 
have two polarization states, + and x . Consider a single Fourier mode with + polarization and as before choose a 
coordinate system with z in the k direction. From Eq. (7.1) in Ref. the contribution of this k-mode to the 
temperature anisotropy is 



T(k, n)+ = ^ J^i 2 3 + ' 6) sin 2 9 cos 20 A+(fc), 



(4.27) 



(note that the choice of the zero of <fi is arbitrary and inconsequential), and A7~ f (fc) is the perturbation to the photon 
brightness induced by this tensor mode after the Polnarev change of variables |56|. For the x polarization state, 
simply replace cos 20 with sin 20 and AT", with A^. Again, we expand this anisotropy pattern in spherical harmonics, 



T,+ 



(k) = i^(2j + l)A+(fc) J dnP 7 (k.n)r ( : m) (n)sin 2 0cos20. 



(4.28) 



Note that for k = z, Pj (cos 9) = [4tt/(2Z + l)] 1 / 2 ^^). Then Eq. flB15| ) can be used to express the integrand in 
Eq. ( 4.28 ) as a sum of products of two spherical harmonics. Orthonormality of spherical harmonics then gives 



l (lm) 



(k) = -(<5 m , 2 +<S m ,_ 2 )^(2j + l) 

3 

, 2Vj(j + 1)0" - l)(j + 2) 



47T 



2j + l 13 



(j + l)(j + 2)(j + 3)(j+4) 
(2j + 1)(2j + 3) 2 (2j+5) 



1/2 



j(j-l)(j-2)(j-3) 



(2j-l) 2 (2j + l)(2j-3) 



(2j-l)(2j + 3) 

The Six's project out only three terms in the sum, and this reduces to 

A+_ 2 (fc) 



1/2 



(4.29) 



x T '+ 



1 



(k) = - Jn{2l + l)(5 m2 + S m - 2 )\ 



'{1 + 2)1 
(l-2)\ 



2A+(fc) 



A/,; +2 (fc) 



(2Z-1)(2Z + 1) (2Z + 3)(2Z — 1) (2Z + 3)(2Z + 1) 



(4.30) 



For the x polarization state, simply replace (5 m2 + <$ m ,-2) with — i[5 m 2 — <$m,-2), and A + with A J 
The contribution of this k mode to C Z T is then 



c ' T ' +(k) = 27TT^ |a ^) (k)|i 



tt (Z + 2)! 
8 (1-2)! 



2A+(fc) 



A Id+2 (k) 



(2Z-1)(2Z + 1) (2Z + 3)(2Z — 1) (2Z + 3)(2Z + 1) 



(4.31) 



and the result for the x polarization state is the same (with the replacement H > x , of course) . If the spectrum of 

+ and x states is the same (which is demanded by statistical isotropy), then the total contribution of tensor modes 
to the temperature anisotropy is 



a 



■T, tensor 
I 



= 2 



d 3 k 



(2tt) 3 1 

1 (Z + 2)! 
8tt (Z - 2)! 



^< + (k) 
k 2 dk 



A|,_ 2 (fc) 



2A+(fc) 



A/,i+ 2 (fc) 



(2Z-1)(2Z + 1) (2Z + 3)(2Z-1) (2Z + 3)(2Z + 1) 



(4.32) 



which agrees with the results of previous calculations [0,0, 



1G 



2. Polarization Moments 



Finally, we calculate the polarization moments from tensor perturbations. The Stokes parameters induced by a 
single tensor Fourier mode with + polarization in the direction n = (9, 0) are J13| 



Q(k, n)+ = -j J^(2j + l)Pj (cos 0)(1 + cos 2 6) cos 20 K+ J (k), 



U(k,h) + = -j-^2(2j + l)P j (cos9)2cos9sin2(f)A^ j (k), 



(4.33) 
(4.34) 



where again Agj(fc) are Legendre coefhcients of the photon polarization brightness for tensor metric perturbations 
with + polarization. Note U is the opposite sign from Ref. [jl3| because the coordinate system there has opposite 
orientation from the one here. The polarization tensor is thus 



pf(k, n) = ^ £(2j + 1) A+ .(*) Mf(n), 



(4.35) 



with 



Af$(fi) = Pj(cos< 



(1 + cos 2 6>) cos 20 -2 cot 6> sin 20 

-2cot6>sin20 -(1 + cos 2 9) esc 2 9 cos 20 



(4.36) 



For x polarization, make the replacements cos 20 — > sin 20 and sin 20 — > — cos 20. 



Now we calculate the multipole moments induced by this single tensor mode. From Eq. (2.17) above, 



(4.37) 



Calculation of this M?Ps :a b is straightforward using Eq. ( AlS ), but the algebra is considerably longer than the scalar 
case. The result is 



M ( a / ):Qb = cos 20 [12(1 - x^P^x) + 8x(l - x 2 )P^(x) + (1 - x 4 )P?(x)] , 



(4.38) 



with x — cos 9. The x pe rturbations give the same result with cos 20 — > sin 20 . Beca use o f this azimuthal dependence, 
the integral in Eq. (4.37) is non-zero only for m — ±2. Using Eqs. (pTJ), (jBq), and (B16), we get 



M ( a / ):afc = cos 20 



(j + 3)(j+4)P^ 2 6j(j + l)P? (j-2)(j-3)P 2 _ 2 



(2j + l)(2j + 3) (2j + 3)(2j-l) (2j + 1)(2j-1) 



(4.39) 



This is just a sum of the three spherical harmonics with m = 2 and the three with m — —2. The integral in Eq. 
(p7|) is 



dnr (im) (fi)M ( a / ):ah (n) = 



/ 7T (j + 2)! 
21 + 1 (l-2)\ 



(o~m2 + 8 m ,~2) 



(l + l)(l + 2)S hj+2 6l(l + l)5 tj 



(2Z — 3)(2Z — 1) (2J + 3)(2J-1) (2Z + 5)(2Z + 3) 

(4.40) 



which results in 



x G '+ 

\lm) 



(k) = -(5 m2 + 6 m ,-2)y/2n(2l + 1) 



(Z + 2)(/ + l)A+,_ 2 (fc) 6Z(/ + 1)A+ (fc) Z(i-l)A+ i+2 (fc) 



(2^ - 1)(2Z + 1) 



(2Z + 3)(2Z — 1) (2/ + 1)(2Z + 3) 



(4.41) 



As before, each k mode and each polarization state contributes to C; in the same way, so integrating over k and 
multiplying by two (for the two polarization states) gives 



a 



G, tensor 
( 



1 

167T 



k 2 dk 



(Z + 2)(I + l)A+,_ a (fc) 6Z(/ + l)A+(fc) 1(1 - 1)A+ l+2 (k) 



(21- 1)(2Z + 1) 



(2Z + 3)(2Z — 1) (2Z + 1)(2Z + 3) 



(4.42) 
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Therefore, (7p> tcni!0r j s due to three Legendre expansion coefficients, Aq,(&), as opposed to an infinite sum as in the 



scalar case. 



The derivation of the multipole moments for the C harmonics is similar. From Eq. (2.18) above, 
«^)( k ) = Y E^' + 1 ) A Q,( fc ) / d{lY (irn^) M^ :ac (h)e c b . 



(4.43) 



This time we get for the "+" gravity wave polarization 



M ( °/ ):ac e c 6 = sin 20 [-8(1 - x z )P' 3 (x) - 2x(l - x*)P$'( 



2\ r>ii 



„2\ u"/ 



sin 2<fi 



2(i + 3) 2 (x) 2(j-2) 
2, • ! ' ' 1 ' 2./ • I ' * 



where the second equality uses the identities (B6) and (B7). For x polarization states, replace sin 20 
The multipole coefficients are then 



* c,+ (k), 



(l + 2)A+(k) + (l-l)A+(k) 



(4.44) 
cos 2(f>. 

(4.45) 



The x perturbations give the same result except for a minus sign between the Kronecker delta's and an overall factor 
of i. Again assuming equal contributions from both + and x tensor perturbations, 



C, c ' tonsor = — I dkk 2 

47T 



(4.46) 



This calculation verifies our qualitative arguments that tensor modes will produce a "C" polarization field. 



Cross Moments 



From Eqs. ( 4.3Cl| ) and (4.41), the non-zero cross moments are 

>W 2A+(fc) 



a 



TG, tensor 



1 



k'dk 



(2/-l)(2i + l) (2/ + 3)(2Z-l) 



8irNi 

(i + 2)(J + l)A+ (fe) , 6i(/ + l)A+(fe) , 1(1 



A|, +2 (fc) 
(2Z + 3)(2Z + 1) 

l)A^ ;+2 (fc) 



(2Z-1)(2Z + 1) 



(2Z + 3)(2Z — 1) (2Z + 1)(2Z + 3) 



(4.47) 



assuming equal contributions from the + and x tensor perturbations. Note the N, 1 in the prefactor comes from the 
temperature coefficient, not from the polarization coefficient. 

Both the temperature and G multipole coefficients for a k tensor mode with + polarization, (k) and a^'^ (k) 
are proportional to (5 m2 + <5 m .- 2 ) (and similarly for x modes with the sum replaced by a difference. On the other 



_ 2 ) (and the same with the 



hand, the corresponding C multipole coefficients, a,^l(k) are proportional to (6 m 2 — 
difference replaced by a sum for x modes). Therefore, after summing over m, the cross moments £7 ; TC > tcnso1 anc ) 
(jGC, tensor ^ ^ enSQr perturbations vanish. This is a consequence of the symmetry arguments mentioned above. 



C. Summary 

In this Section we have calculated the CMB temperature and polarization multipole moments for both scalar and ten- 
sor perturbations. For theories in which the perturbations have a Gaussian distribution (such as inflationary models), 
this set of multipole moments fully specifies the statistical properties of the combined temperature/polarization map. 
In virtually all theories considered, scalar and tensor perturbations are statistically independent, so their contributions 
to the CMB temperature anisotropy and polarization add in quadrature. Even if they are not statistically independent, 
angular orthogonality of the different modes (i.e., scalar, vector, and tensor modes in the z direction induce non-zero 
multipole coefficients only for m = 0, m = ±, and m — ±2, respectively) guarantees that the contributions of scalar, 
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vector, and tensor perturbations to the multipole moments will add in quadrature. Therefore, for Gaussian theories, 



the map will be fully described by Cj- = C l 
[cf., Eqs. ( E~16|) and (jy|)], 



T.scala 



a 



T, tensor 



a 



C, tensor 



[cf., Eqs. <K§ and (4.32)] 



[cf., Eqs. (4.46)1, and Cj G = C* ( 



G, scalar 



GG /^G, scalar . ^G, tensor 



a 



TG, tensor 



[cf., Eqs. (|4~26| ) and 



(4.47)]. For non-Gaussian theories, n-point correlation functions with n > 2 may be nonzero. 

Eq. (4.4) for the temperature moments from scalar perturbations is written as an integral of a square of a single 
An. However, all of the other moments are written as squares of a sum of several Aji's and/or Aq;'s. This is because 
a spherical-harmonic decomposition is not natural for a Stokes parameter Q, which is a tensor component, nor for the 
brightness perturbation from tensor modes. The calculation can be reformulated using a tensor-harmonic expansion 
for these quantities, which provide a natural basis |BU. 



D. Line-of-Sight Approach 



A very efficient and accurate algorithm for computation of multipole moments has recently been given by Seljak 
and Zaldarriaga |3^] and applied to the polarization multipole moments [^lj (hereafter, SZ). Although significantly 
different in appearance, their spin-harmonic formalism is equivalent to our tensor-harmonic formalism. Here, we 
briefly compare our results with theirs. Although the formalisms differ and the calculations are somewhat lengthy, 
we find that the results agree, which gives us confidence in both sets of results. Furthermore, by identifying the 
moments in their paper with those in ours, their numerical code (which has been made publicly available) can be used 
to compute the multipole moments presented in this paper. 

Consider, for example, the G polarization moments from scalar perturbations. According to Eqs. (12) and (14) 
in S Z, th e polarization pattern induced on the sky by a scalar k mode in the z direction (i.e., their version of our 



Eq. tyL§) is 

O /•TO 

Q(n) = |(1 - cos 2 9) J dr e lkTCOs9 5 (r)H(fc, r)£(k), (4.48) 

and U(n) = 0, where r is the conformal time, g(r) is a visibility function, II is a combination of intensity and 
polarization perturbations, and £(k) is an amplitude for the scalar mode (see SZ for more details). Therefore, the 
polarization tensor induced by this scalar mode is 

P ab (h) = ^(k) jH dT g (r)U(k,r)e^ cose ( J _ Q ) . (4.49) 

Using the rules of covariant differentiation, we find that 

cos 6 d 



P ab :ab (n) = ^(k) £° dr 5 (r)H(fc, r) (JL + 3 



sin(9 d6 

•2 



2 e l 



= drg(rMk,r)^^ [(1 - cos 2 0) 2 e ife " os9 ] , (4.50) 

which agrees with their A^p/ 2 given in their Eq. (15). The C l ' scaal moments are obtained by plugging this into 



Eq. (2.17), squaring, summing over m and then integrating over k. We then find that our results agree with theirs 
(realizing that our (27r) 3 Fourier conventions differ from theirs) if we identify Cp = Cei/2, where Cei are their 
electric-type moments. We have further checked that our temperature moments agree with theirs (with no factor 
of two), and that our C moments are twice their B moments: Cf — Cbi/2. Our Cj G are equal to their Cci/V^- 
Although we do not present it here, we have checked their tensor-mode calculations as well. Alth ough signifi cantl y 
more involved, we still find agreement; for example, compare their £(x) and B(x) with our Eqs. ( 1.38 ) and ( |l.44| ). 



The identification of our polarization moments with theirs is also consistent with our Wp 1 in Eq. ( 3.18 ) being half 
theirs. 



V. TWO-POINT CORRELATION FUNCTIONS 

In this Section, we relate the multipole moments C; T , Cp, Cf, and Cj G to two-point correlation functions of 
temperature and Stokes parameters. This will make contact with previous work on the subject. We also derive 
fiat-sky limits useful for analyzing maps of small sky patches. 
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A. Correlations Between Temperature and Stokes Parameters 



The linear-polarization state at any given point is specified completely by the Stokes parameter Q and U, but these 
quantities depend on the coordinate system which one chooses. On the other hand, we know that Q and U transform 
as the components of a STF 2x2 tensor, so given Q and U in some coordinate system, we can always determine Q' 
and U' in any other coordinate system. 

The Universe is assumed to be statistically isotropic, so it is possible to construct two-point correlation functions 
which depend only on the angular separation between the two points. But simply correlating Q and U in a particular 
coordinate system gives correlation functions which depend on the positions of the points being correlated as well the 
angular separation. This is what has been done in previous published work. 

A coordinate-independent set of correlation functions can be expressed in terms of the ones which have appeared in 
the previous literature. The prescription is simply to define correlation functions of Stokes parameters with respect to 
axes which are parallel and perpendicular to the great arc (or geodesic) connecting the two points being correlated. 
So Q r is the difference in intensities in two linear-polarization states parallel and perpendicular to the great arc 
connecting the two points, and U r is the difference in two linear-polarization states which lie 45° away from the 
parallel and perpendicular. The three quantities T, Q ri and U r have six correlation functions between them: (TT), 
(U r U r ), (Q r Qr), (QrT), (Q r U r ), and (U r T). However, only four can be nonzero. Although Q r and T are invariant 
under reflection along the great arc connecting the two points being correlated, U r changes sign. Therefore, the 
expectation values (Q r U r ) and (U r T) must be zero from statistical isotropy. This is as expected: four nonzero sets of 
moments Cj , C\ • , C\ p , and Cj describe the map. Correspondingly, four nonzero correlation functions provide an 
equivalent statistical description. 

We begin with the familiar (TT) correlation function, 

C^) = /^l^)\ . (5.1) 



To To 



ni -ri2— cos 6 



The correlation function depends only on the angular separation of the two points, so in calculating it, we may 
choose one point to be at the north pole, (0,0), and the other to be on the <f> = longitude at a distance 8 from 



the north pole, (6, 0). Then expand T(n) in terms of spherical harmonics as in Eq. (2.8) and note that Y(/ m )(0, 0) 

T(O,O)T(0,O)' 



^/(2Z + 1)/(4tt). So 

C T (8) 



To T 

E ( a Ilm) a Vrn>) Y *l m ) («, 0)l"(J'm') (*> °) 



Iml' 



ItnV 



Air 
I 



where we have used Eq. (3.1) to go from the second to the third line. This recovers the well-known result for the 
temperature autocorrelation function. 

The derivation of the polarization correlation functions will proceed analogously (and is similar to the case for 
weak-lensing correlation functions p8| ), and requires Wn m \ and Xn m ) at the north pole. Using the asymptotic 
relations 

/i\(m+H)/2 (i , i m m 

p !"^~ 2 H H ! I- -Mil ( "^°» (5 - 3) 



Pi(cos6) - 1 - h(l + 1)6 2 , 6^0 (5.4) 

it is straightforward to show that for \m\ > 2, Xn^ and Yn m ) are both asymptotic to 6>l m '~ 2 as 9 — > 0, so they are 
nonzero at 8 = only for \m\ — 2; for m — and m = 1, they are asymptotic to 2 and 8. After a little algebra, 
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W (lm) (0,0) 



1 21 + 1(1 + 2)1 

2 V 4tt (7-2)! 



(<5m2 + <5m,-2), 



and 



i 2/ + l([+2)! 
A (im )(0,0) - 2 V (? - 2)! ^ ~ dm '- 2 >- 



Now consider the (QQ) correlation function 

C Q {9) - 



Q r {ni) Qr(n 2 ) 



(5.5) 



(5.6) 



(5.7) 



ni -ri2— cos t 



where, once again, the Stokes parameters Q r are defined with respect to axes parallel and perpendicular to the great 
arc connecting fix and ri2- As in the temperature case, choose one point to be at the north pole and another a 
distance 9 away along the = longitude. This choice has the added advantage that the great arc connecting these 
two points is along the 9 direction, so we can use the Q defined in the (9, 4>) coordinate system, Eq. (3.2). However, 
in this coordinate system, the definition of Q at the north pole is, strictly speaking, ambiguous. Therefore, we always 
consider a point on the (f> — lon gitu de which is infinitesimally close to the north pole; in other words, Q(0,0) really 
means lime ^, Q (9,0). Using Eq. (3J3) for the Stokes parameters, the (QQ) correlation function is then 



C Q (9) = 



/Q(O,O)Q(0,O) 



To 



E N l N i'{l*?im)W ilm )(0,0) - a^„ l) X ( , m) (O,O)][a^4, ) lT ( * i , m , ) (0,O) - af^Xfr^B, 0)}) 

Iml'm' 

E \f^^ N ^( W m + W (i,-2)) + iC?{ X h) ~ X (*,-2))L (5-8) 



where we have used Eqs. (3.1), (5.6), and (5.5), and (.Q?im) a ?im) ) = ^- This can be simplified using X^ 



and Wn m \ = Wn _ m N and the definitions in Eqs. ( 2.24 ) and ( ^.25| ), giving 



L (im) 



-X 



C Q (9) = E -^ r N ?[C?G+ 2) (cos9) + CfGfo (cos 0)] 
i 

For the (UU) correlation function, the derivation is similar, giving 

= -Y,^ N ?[CFGl l2) (cos9)+Cl i G m (cos9)}. 
i 

For the (TQ) cross-correlation function, 

TQ = / T(nt) Q r (n 2 ) 
\ T T 



(2,-m) 



(5.9) 



(5.10) 



11 1 II o —COS f 



E^^C-P^eos, 



(5.11) 



Eqs. (5.9), ( 5.10D , and (5.11) are exact (i.e., there is no small-angle approximation) expressions for the polarization 
correlation functions. 



B. Multipole Moments from Correlation Functions 



Above, we derived expressions for correlation functions in terms of multipole moments, and now we perform the 
inverse transform and express the multipole moments in terms of correlation functions. Begin with the temperature 
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autocorrelation function: multiply both sides of Eq. ( p.2[ ) by the Legendre polynomial Pi'(cos9), integrate over cos 9, 
and use the orthogonality of Legendre polynomials to obtain 



Cj = 2tt / dO sin 9 Pi (cos 9) C T (9) . 
Jo 



(5.12) 



Similarly, for the polarization-temp erat ure moments, multiply both sides of Eq. ( |5.11 ) by the associated Legendre 
function P 2 , integrate, and use Eq. (B3) to obtain 



Cf G = irNi f d9 sin BP? (cos 0) C* TQ (9) . 
Jo 



(5.13) 



The derivation of the polarization moments from the polarization autocorrelation functions is similar. Orthonor- 
mality of the tensor harmonics implies that 



J {W ( * lm) (n)W (l , ml) (h) + X* (lm) (h)X (l , ml) (h)} dn = -^5 w 5 mm >, 



and it can also be shown that 



-*(* m )( fi )W(i'm')(fi) + W *lm)(n) X (l>m>)(n)} d& = 0. 



(5.14) 



(5.15) 



To do so, note that the 4> dependence of W(i m \ and Xn m -\ is just e tm ^, which means that the integral is immediately 

r (lm) 



zero for m^m'. For m — vn! , the integral over cos9 vanishes using the explicit forms of Gf, |28|. From Eqs. (|5.9|) 



and ( 15.101 ), 



[C Q {6) + C u {6)]e 21 * =Y, 



21 + 1 
2tt 



iV,(Cf + Cp)[W {l2) (9, 0) + iX (I2) (0, 0)] 



(5.16) 



[C <3 (6>)-C c/ (6')]e 2i<?i =^ 



2tt 



^(Cf- - Cft^CM) - iX {m {9,4>)}. 



(5.17) 



Then multiply both sides of the first equation by W* l2 ^ — iX? l2 s and the second by W* l2) + iX* l2) , integrate over all 



directions n, and apply Eqs. ( 5.14 ) and ( 5.15 ), giving 



(12) ^ l ^(l2)' 



c? + cf = 



2vr Ni 

2T+Y2 



dh[CQ{9) + C u (9)} e 2 ^ [^ 2) (n) - iXfafr)], 



(5.18) 



C? Cf 



2tt Ni 



dn[CQ(9) - C u {9)] e 2 ^ [W( l2) (n) + iXfa{&)]. 



21 + 1 2 

Upon summing and differencing and carrying out the integration over cj>, we obtain 

CP = -TriV 2 / d6»sin6»[C Q (6i)G+ JCOS6*) +C t/ (6»)G" 2l (cos6')], 



(5.19) 



(5.20) 



TriV 2 / d9 sin 9[C u (9)G+ 2) (cos 9) +C Q (9)G~ 2) (cos 9)} 



(12) 



(5.21) 



which are the desired relations giving the polarization multipolc moments i n ter ms of the polarization autocorrelation 
functions. Given some measured correlation functions, evaluation of Eq. (5.21) for any / will probe the existence of 
non-scalar modes. 
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C. Correlation Functions in the Small-Angle Limit 



In order to make contact with previous work and to present estimates useful for measurements on a small patch 
of the sky, we now derive the small-angle limit of the above expressions which give correlation functions in terms 
of multipole moments (and vice versa). The correlation functions given in previous work were of Stokes parameters 
measured in a fixed coordinate basis, whereas ours are of Stokes parameters measured with respect to the great 
arc connecting the two points being correlated. However, our results can be compared with previous results by 
taking tf> = in their expressions. Although the expressions for correlation functions in the small-angle limit given 
in Refs. |T3"l|l4|| are quite complicated, when the small-angle limit is taken consistently in all steps, the expressions 
simplify greatly, as emphasized in Ref. [p9|l ( and resemble correlation functions for ellipticities of galaxies due to weak 
lensing from large-scale inhomogeneities ||28|j40| ) . 

Once again we begin with the temperature moments. A useful asymptotic relation is 



Pi (cos 9) ~ J (fl), 



where J m (s) is the Bessel function of order m. Substituting into Eq. (|5 
taking the limit I 3> 1 gives 



s = (21 + l)sin(0/2) -> 0, (5.22) 
[), approximating the sum by an integral and 



(5.23) 



C T (9)~—j ldlJ (W)C? 



for 9 <C 1. For the temperature-polarization cross-correlation function, we note that 

if (cos 9) ~ 4/ 4 J 2 (s), s^O, 

which gives 

2 3/2 



C T ®(9) 



PdlCj G J 2 (W), 



for 9 <C 1, from Eq. ( pUl]) . 

For the polarization autocorrelation functions, note that 



(5.24) 



(5.25) 



(COS 9) 



l [J (s)± Ms)], 



From Eq. ( |5.9| ) we obtain 



and from Eq. (5.10) we obtain 



ldl[(Cf 



ldl[(Cf 



C?)M19) 



C?)U19) 



ic? 



c?)MW)], 



Cp)J 4 (l9)}, 



for 9 < 1. 



(5.26) 



(5.27) 



(5.28) 



Eqs. fl5.23| ), ( |5.25| ), ( |5.27| ), and (|5.28[ ) agree with the forms in Eq. (19) in Ref. |3g] for 4> = and Cf = 0. (Also 
recall that (Q r U r ) = (TU r ) = 0.) If Cp = 0, then (Q r Q r + U r U r ) (which depends on an integral over Jo (10)) 
and (QrQr — U r U r ) (which depends on an integral over J ±(10)) depend on the same set of moments Cp and are 
therefore not independent. However, if Cp 7^ 0, then these correlation functions will depend on two independent sets 
of moments. 

We can also derive expressions for the multipole moments for / 3> 1 in terms of a correlation function measured at 
small angular separations. For example, using Eq. (5.22) to approximate Eq. ( |5.12 ) for / 3> 1 gives 



C 



2tt / J o (W)C T (0)0d0 



and using Eq. ( [3.24 ) to approximate Eq. ( 5.1 3| ) gives 



TG 



Ay/2wl 2 / J 2 (W)C TQ (0)0d0. 



(5.29) 



(5.30) 
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Using Eq. ( 5.26| ), we can approximate Eq. ( |5.20| ) by 



Cf ~ -| / (C Q (0)[J o (/0) + J 4 (W)] + C^MJo^) - J 4 (W)]) Ode, (5.31) 

Cf ~ -| / (C u (9)[J {W) + J 4 (/0)] + C Q {9)[MW) - J±(W)]) OdO. (5.32) 

If any nonzero Cp is found in this way with correlation functions measured on a small patch of the sky, it is an 
indication of vector or tensor modes. 

VI. SUMMARY AND DISCUSSION 

This paper provides a detailed and complete formalism for characterizing polarization fluctuations in a full-sky map 
of the cosmic microwave background. We give explicit forms for tensor spherical harmonics, in which the polarization 
can be expanded in direct analogy to the expansion of temperature perturbations in the usual spherical harmonics. 
The tensor harmonics are numerically just as easy to evaluate as spherical harmonics, so polarization map simulation 
and analysis will be no more cumbersome than in the temperature case. 

The most important physics results presented here are that of the six sets of multipole coefficients describing the 
correlations in a temperature/polarization map, two must be zero if the Universe is parity-invariant, and a third 
vanishes for scalar metric perturbations. The moments Cp , which are non-zero only for vector and tensor metric 
perturbations, are in principle an unambiguous probe of primordial gravity waves and vorticity JlS| , |20) ] . A cosmological 
contribution to the moments C ; TC or Cp would demonstrate a remarkable handedness to the primordial perturbation 
spectrum. A much more likely and practical use of these moments is to monitor foreground microwave emission. We 
also note that CMB polarization may give useful information on primordial magnetic fields Jl7],^l| and galaxy cluster 
magnetic fields p2[ . Measurement of polarization in the Sunyaev-Zeldovich effect can be used to measure cluster 
transverse velocities jl3| and/or the CMB quadrupole moment incident on the cluster ji^] . 

Most current microwave background codes calculate the Legendre coefficients of the radiation brightness in Fourier 
space [ p5[ . We have derived exact expressions for all of the multipole moments in terms of these brightness coefficients. 
For tensor metric perturbations, the expressions are particularly simple and trivial to implement numerically. The 
result for scalar perturbations is somewhat more complex, involving an infinite sum over the brightness moments. 
However, the contribution of the sum to the final expression for the multipole moments is only significant for the 
lowest moments, so the overall cost of the computation should only increase slightly. The formulas for the multipole 
moments derived in this paper should allow for relatively simple conversion of existing CMB codes. 

Of course, a cosmological signal will have to be distinguished from foreground contamination. Synchrotron emission 
from our galaxy is highly polarized and extragalactic radio sources may also contribute significantly ^6|. The 
amplitude of these foreground polarization sources is unknown at the present time. Since both likely foregrounds have 
a spectral dependence substantially different from the blackbody CMB spectrum, the usual techniques for subtracting 
foregrounds from temperature maps should also work for polarization WA . 

Of course, simply attaining the necessary sensitivity to make any polarization detection will be a great experimen- 
tal accomplishment. The MAP satellite, currently being constructed, will have the sensitivity to make a statistical 
detection of polarization. The COBRAS/SAMBA satellite, now in the planning stage, should be capable of seeing 
polarization on a pixel-by-pixel basis if it is configured to measure polarization. At this time, it is undecided whether 
COBRAS/SAMBA, which uses incoherent bolometer detectors in its most sensitive frequency channels, will sacrifice 
some temperature sensitivity to make polarized measurements. But optimistically, within a decade we may have in 
hand detailed temperature/polarization maps of the cosmic microwave background. How much cosmological informa- 
tion can be extracted from such maps is currently under study |m| . The formalism presented in this paper provides 
a basis for addressing such questions. 



and Eq. (5.21) by 



ACKNOWLEDGMENTS 



We would like to thank Jerry Jungman, Lloyd Knox, and Uros Seljak for helpful conversations, and David Spergel 
and Gary Hinshaw for stimulating questions. This work was supported by D.O.E. contract DEFG02-92-ER 40699, 
NASA NAG5-3091, and the Alfred P. Sloan Foundation at Columbia, NASA AST94-19400 at FNAL, and the Harvard 
Society of Fellows. M.K. acknowledges the hospitality of the NASA/Fermilab Astrophysics Center and the CERN 
Theory Group. 



24 



APPENDIX A: DIFFERENTIAL GEOMETRY ON THE SPHERE 



This appendix collects results from differential geometry, with particular application to the manifold S 2 (the 2- 
sphere) , which are needed in definitions of and calculations with the tensor spherical harmonics on the celestial sphere. 
We use the notation f, a = df/dx a to indicate a regular partial derivative and f- a = V / for a covariant derivative. We 
use the colon, ":" , rather than the more traditional semi-colon, ";" , to distinguish derivatives on S 2 from 4-dimensional 
derivatives in general relativity. All of our tensors are defined with respect to a coordinate basis. Note the the metric 
tensor g a b commutes with covariant differentiation: g a b.c = 0. The determinant of the metric is denoted by g = \\g a b\\- 
Covariant derivatives of scalars, vectors, and (rank-2) tensors are 

S.. a = S, a V a :h = V\ b + V c T a bc T ab :c = T a \ c + T db T a cd + T ad T b cd (Al) 

where the T^ c are Christoffcl symbols defined by 

Kc = \ 9 ad (9db,c + 9dc.b - 9bc.d)- (A2) 

We have used the formula for the 4th derivatives of a scalar function 

S' ab a b = V 2 V 2 S' + R dh S : db + 2 R' d S.d W 2 S = S' a a R a b = R c ac b R=R a a (A3) 

where Rabcd, Rab, and R are the Riemann tensor, the Ricci tensor, and Ricci scalar, respectively. We have also used 
the formula of integration by parts 

j> d 2 h^X ab Y ab = -l d 2 h^X ab :a Y.. b = I d 2 h^X ab .. ab Y, (A4) 

where § indicates integration over a closed manifold with no boundary, and n is shorthand for dx x dx 2 . 

On a two dimensional manifold the Levi-Civita symbol is a traceless anti-symmetric rank-2 tensor given by 

eab = \/9^_° 1 J) (A5) 

which has the following simple properties: 

£cae C 6 = g a b = —£ac£ C b, tabbed = 9ac9bd ~ 9ad9bc, ^ab:c = 0. (A6) 

A rank-2 tensor is trace- free (or traceless) iff g ab T a b = 0, and a rank-2 tensor is symmetric iff T a b = Tb a - In two 
dimensions the latter requirement is equivalent to e ab T a b — 0. Linear polarization is described by a symmetric 
trace-free rank-2 tensor, or STF for short. 

The geometry of a two dimensional manifold is particularly simple because it is determined solely by the Ricci 
scalar R. Some handy identities are 

Rabcd = ^ R tabbed Rab = ^ R9ab ^ ' Rabcd = R^cd ^ Rabcd = ^R^bd- (A7) 

Another useful identity is 

M ab N cd e ac e bd = ~M ab N ab if g ab M ab = g ab N ab = e ab M ab = e ab N ab = 0, (A8) 
i.e. for STF tensors only. 

In this paper we are only interested in the geometry of the unit sphere, S 2 . Its geometry is nearly as simple as 
can be, since R = 2. We exclusively use the explicit coordinate system defined by spherical polar coordinates, (9, 4>), 
where 9 is the polar angle from a particular point on the sphere, and <j> labels the angle on circles which are centered 
on this same point. In these coordinates the metric is 

9ab = ( J gin °2 B ^ , g=\ | gab ! I = sin 2 9. (A9) 

while the antisymmetric tensor is 
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/ sin#\ a ( sin#\ ,., n , 

e "»= V - S in0 J' eb ={-c S c6 J" ^ A1 °) 

From this metric the Christoffel symbols follow as 

r^ = -sm0co S r^ = rJ e = cote r e ee = v e H = r% = vt e = v% = o. (An) 

The explicit components of the 2nd derivatives of a scalar function are 



Y :ab 


= Y ab 


Ym 


= Y es 


Y.. g<t> 


= Y M 


Y-ss 


= Y.64 



cot e y„ 



(A12) 



while explicit expressions for the some second derivatives of a symmetric rank-2 tensor are 

M ab :ab = M ee . gg + 2M°+ M + M** M - sin0 cos fi 

+ 2cot0M 8 % + 4 cot (9M e % + (1 - 3 cos 2 0)M*+ - M ee (A13) 

and 

M ab :ac e c b = sm6 (M e +,oo + M++ - csc9(M ee M + M* 6 ' M ) - cot 9 esc 9M ee ^ 

+5 cos 9M e ^ fi + 3 cos 9M ^ + 3 (cos 6 cot 9 - sin 9) M e<i > , (A14) 

which again only apply if e ab M ab = 0. 

APPENDIX B: SPHERICAL HARMONIC AND LEGENDRE FUNCTION IDENTITIES 

In this Appendix, we list some identities involving spherical harmonics and Legendre functions which we have used 
in our calculations. 

The associated Legendre functions are defined by 

pr{x) = (-ira - x *r/ 2 — Pl ( X ), m > o, (bi) 

p- m (x) = (-i r f—!^p r . (B2) 
(I + my. 



Normalization: 



Legendre Polynomials: 



1 9 (7 + roV 

dxP{ n {x)Pr{x) = 6 U > , , , I (B3) 



i 



21 + 1 (l-m)\ 



Pl{x) = P^{x). (B4) 
The associated Legendre functions satisfy the following recursion relations (see Eqs. 8.733 in Ref. fljsfl): 

P/" +2 (x) + 2(m + 1) X P™ +1 M + {l~m){l + m,+ l)P ; m (x) = 0, (B5) 
VI — x 2 



(21 + lJvWPr'W = - iftiO*)- (B6) 

(2/ + l)xP/" = (2 - m + l)P^i + + m)PJlx, (B7) 



2G 



dp[" 

dx 



(m+l)xPr-(l-m + l)P^, 



(B8) 



(1 - x 2 )-^r = ~^Pi l + (I + m)Pl\. 



(B9) 



Using orthonormality and completeness of Legendre polynomials, an associated Legendre function of m = 2 can be 
written in terms of Legendre polynomials using the integral |49[ 



dxPi{x)Pf{x) 



0. 



I > j or j + I odd; 



-2j(j-l)/(2j + l), j = l; 

4, I < j and j + I even 



With the above conventions for Legendre functions, the spherical harmonics are given by 

WM) = J 2l l 1 n~ n l ) !', Pr(cose)e 



Air (l + m)r l 

The convolution of a Legendre polynomial with a spherical harmonic is 



im0 



dhPjCk- n) y (Zm) (n) 



47T 



2/ 



(BIO) 



(Bll) 



(B12) 



This can be obtained by expressing the Legendre polynomial in terms of spherical harmonics with the spherical 
harmonic addition theorem and then using the orthonormality of spherical harmonics. 

Angular-momentum lowering and raising operators can be used to derive the following recursion relations for 
spherical harmonics: 



^sm0Y {lm) = 



(I + m + l)(l + m + 2) 



(2Z + 1)(2Z + 3) 



1/2 



*(i+l,m+l) + 



(I — m)(l — m — 1) 



(2Z-1)(2Z + 1) 



1/2 



Y(l-1 



m+1) 



(B13) 



(h, 



(I -m + l)(l -m + 2) 



1/2 



Y(i+i 



(2Z + 1)(2Z + 3) 
which can be iterated and evaluated for m — to obtain 



(Z+l,m-l) 



(I + m)(l + m - 1) 
(2Z-1)(2Z + 1) 



1/2 



Y(i-i 



(l—l,m—l)> 



2 1 

cos 2cj) sin 6Y( l0 ) = - 



(l + l)(l + 2)(l + 3)(l + A) 



n 1/2 



(2Z + 1)(2Z + 3) 2 (2/ + 5) 



(X(i+2,2) + Y( l+2 .-2)) 



y/Ki -i)(/ + m + 2) 

- 2 — - (Y(L2) + Y (l-2)) 



+ 



(2Z-l)(2Z + 3) 
l(l-l)(l-2)(l-3) 11/2 



(21 -3)(2Z- 1) 2 (2/ + 1) 



(^(1-2,2) + ^(J-2,-2)) 



(B14) 



(B15) 



For the replacement cos 20 — ► sin 20 make the replacements (1/2) — * (l/2i) in the prefactor and Wx.a) + ^(^,-2) 
(X(a;,2) - ^(x -2))- 



Looking at the 8 dependence of Eq. ( B15 ), we get 
(l-x 2 )P 3 (x) 



p2 



2Pf 



pi 



(2j + l)(2j + 3) (2j-l)(2j + 3) (2j + l)(2j-l) 



(B16) 
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